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FOREWORD 


The work described herein was done by The Boeing Company, 
Wichita Division, under NASA contract NAS 3-14321 with Mr. 
H. Bloomer, V/STOL and Noise Division, NASA -Lewis 
Research Center, as Project Manager. 


ABSTRACT 


An acoustic radiation analysis has bean developed to predict the far-field characteristics of fan noise 
radiated from an acoustically lined unflanged duct. This analysis is comprised of three modular 
digital computer programs which together provide a capability of accounting for the impedance 
mismatch at the duct exit plane. This report discusses the Directivity Index Program whose 
relationship with the other two modular reports of the analysis is illustrated on the following page. 

Admissible duct configurations include circular or annular, with or without an extended 
centerbody. This variation in duct configurations provides a capability of modeling inlet and fan 
duct noise radiation. 



DEVELOPMENT OF ACOUSTIC RADIATION ANALYSIS 
OF TURBOFAN NOISE FROM LINED, UNFLANGED DUCTS 

OVERALL REPORT ORGANIZATION 
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1.0 


INTRODUCTION 


1.1 Background 

This report discusses the Directivity Index Program that provides the radiation portion 
of the package of modular computer programs which together form a first generation 
capability for the analytical prediction of far field noise being radiated from an 
unflanged inlet or exhaust duct. 

The development of this program had its origin in the need for a method to determine the 
directional distribution of acoustic energy in the far field which is radiated into free 
space from the fluctuating pressure field appearing at the open face of an unflanged 
duct. 

This digital computer program utilizes the output from the Acoustic Source Distribution 
Program, Reference 1, to generate die far-fietd radiation patterns for a duct with no 
open-end correction. The impedance mismatch at the duct face can be accounted for by 
the Duct Termination Impedance Program, Reference 2, to determine the correct 
pressure distribution on the open duct face. 


1.2 Technical Approach 

This technical development addresses the far-field radiation pattern that is generated 
from fan or compressor noise exiting from an inlet or exhaust duct. Unflanged ducts 
with both annular and circular cross-sections can be analyzed by the program. The 
option of extending the centerbody beyond the duct face permits the program to 
analyze the three duct configurations shown in Figure 1. Any combination of these 
configurations may be superimposed to model a specific engine geometry. For example, 
by representing an engine exhaust duct with Configuration B and superimposing it upon 
Configuration A which represents the inlet, the far-field radiation patterns are calculated 
for the Inlet and exhaust separately and then combined. 

The approach presented is an extension of the case of sound radiating from a duct 
enclosed by an infinite baffle to the unflanged case which accounts for the diffraction of 
sound around the outside lip of the duct. The phased contributions of the outer and 
extended inner duct surfaces to the total far-field pressure must be determined and then 
combined with the radiation of acoustic energy from the duct face. 

The far-field radiation characteristics of the duct configuration are determined from the 
spatial distributions of acoustic sources over the duct walls and exit plane. The source 
disttf.jution which satisfies the imposed boundary conditions is determined by the 
Acoustic Source Distribution Program, Reference 1. These boundary conditions require 
that the acoustic perturbation velocity normal to the outer duct walls must vanish. In 
addition, the normal velocity distribution over the duct exit plane must equal the 
velocity distribution associated with either a single mode or a combination of modes of 
discrete frequency. 
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Given the source distribution over the duct, the far-fieid directivity is found by 
numerical integration of the Helmholtz integrals for the velocity potential (under the 
far-field assumptions) over the duct walls and exit plane. Far-field radiation patterns can 
be determined separately for the inlet and exhaust ducts and then combined. 

Several assumptions have been made in this development. These concern the transport 
characteristics of the model and the decomposition of the acoustic sound energy into the 
acoustic modes present in the duct. The acoustic sound is assumed to be radiated 
through an inviscid perfect gas which has no mean flow. This simplification will not 
significantly affect the results for noise propagation through a locally subsonic flow field 
of an inlet. Radiation of sound through an exhaust or jet flow, a case where the 
refraction of sound due to the velocity gradients is important, will not be represented 
properly by this model. 

The velocity distribution is assumed to have a Cos (mf) ) type of angular dependence. 
This simplification is justified by the angular characteristic functions determined from 
the solution of the governing wave equation for sound propagation through annular 
ducts. This assumption does not affect the ability to couple different modes in the 
far-field program, but does not restrict the Duct Termination Impedance Program, 
Reference 2, to the extent that for a velocity distribution whose angular dependence is 
specified by a combination of different m angular modes, a separate analysis must be 
made for each set of angular modes. 



2.0 


MATHEMATICAL DEVELOPMENT 


The far-field radiation characteristics of an unflanged duct configuration are formulated 
in terms of the acoustic velocity potential. This formulation provides a mathematical 
model with a straight forward method of satisfying the velocity boundary conditions on 
the duct walls and exit plane. Refraction of sound from the duct interior face to the 
far-field is inherently accounted for in the phase relationship of the potential function. 

Since the velocity potential is governed by a linear set of equations, superposition may 
be used to determine the far-fietd potential 0 . The far-field potential for the math 
model shown in Figure 2, is determined by the superposition of three 
potentials, $> 0 , 0 ( , <p 2 , which represent the contribution of the duct face, outer duct 
wall and extended inner duct wall, respectively. This technique is expressed 
mathematically by 


4> ( r f Q r x,i ) * 


(1) 


The component velocity potentials are determined by a spatial distribution of monopole 
acoustic sources on each reflecting surface. For example, the velocity potential due to a 
distribution of harmonically vibrating sources on the face of the duct is 


<k = e 


ff/ 

•'o Ja 


\( S .) e iwR ° /c 


cos(m0 o ) r o d0 o dr Q 


(2) 


The form of this integral equation is derived from the following assumptions: 

• Acoustic velocity at the point P may be determined from a monopole 

acoustic source, ^ which represents an outward 

propagating wave and is a solution of wave equation in spherical coordinates. 

• The angular dependence of the potential may be described in terms of the 
characteristic function Cos (m 6 ) for a given angular mode number m. 

■ The radial distribution of sources is given by the spatial distribution 
function 0 o (r) which is determined by the Acoustic Source Distribution 
Program, Reference 1. 

Similarly, the velocity potential terms for the outer and inner duct wall surfaces, 
respectively are 

r® f, ,p a 

= 6* W I f cos(m0,) d0 b dz 

’ L i R ' 131 
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( 4 ) 



e 


cos(m0 2 ) df? 2 a dz 2 


The classical method of evaluating these integrals is to utilize the far-fietd 
approximations, that is, the distance between the source and the far-field point is much 
greater than the argument of the source strength. 

s » r 0 

The separation distance R 0 in spherical coordinates 
Ro = s 2 -2sr 0 Sin ip Cos (Q - d Q ) 
may be approximated by 

R 0 == s*r Q sin 4> Cos { 6 - 0 Q ) 

In the exponential term and 

R 0 ^ s 

in the reciprocal term. 

The separation distances for the outer and extended inner duct surfaces are 
approximated by expanding the radical expression for R-j and R 2 with the truncated 
binominal expansion. 


R, « V s 2 -2SZ, Cos ^ - bSinlAcos(0 - 0,) 

1 - 1 1 
Rt ~ s 2 -2sz, Cos tf'+z, ~ u, 

R a ~ s 2 - 2 sz a Cos ^ +z| ‘ - aSinl/'Cos (Q - $ 2 ) 

111 
R 2 ■yj a z -2sz a Cos +Zj u 2 
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Utilizing the far-field approximations in the integral Equations (2) through (4) the 
far-field potential of Equation (1) becomes 


4> 


= e 


iu s/c r 

— / / 

•^a Jn 


b ^ 0 Sin ^Cos (0 - 0 O ) 


_ » £J -27T 

+ Afo )e _ _ dzJ e 
J u > 


Cos(m0 o )r o d 


^bSin^Cos(0 - 0,) (5) 

Cos (m 0! )bd0 T 


-CO . ■ 

J <t>^e Lc 2 dz J e 

Jq U 2 Jo 
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Cos (m0 2 )ad02 


A transformation of the trigonometric functions and the definition of the Bessel 
function of the first kind: 



--m C £ z Cos 0 

IT f 
-i r 


Cos (m 0 )d0 


(6) 


permits the 0 integration to be made in closed form. 

Utilizing Equation (6) and collecting terms, the potential function in Equation (1) can 
be written in terms of the directivity functions. 

LlOS/C r- -| 

« 27T [~L) m Cosm0^-^ — |G 0 ( f /') + G,{^)+G 2 C^)J (7) 

The directivity functions which indicate the contribution of the duct face, outer and 
extended inner duct walls to the directional distribution of energy in the far-field are 
defined as 



G 0 {^) = / AM Jm Hr r ° sin ^ ) r o flr o 

Ja 


G, ($) " sb J m (-~-bSin 
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(9) 
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( 10 ) 




3.0 


METHOD OF SOLUTION 


The constraint that this program must be compatible with the Acoustic Source 
Distribution Program of Reference 1, required that the "Box Method" technique be used 
to solve Equations (8) through {10}. This mathematical technique assumes that the 
variation of the source distribution over a small portion or box of the reflecting surface 
is constant, an assumption which allows the integrations to be evaluated by a 
summation. To facilitate the numerical solution of these equations, the duct exterior 
face is divided into a sequence of M annular rings, and the duct walls into N and N2 
cylindrical segments as shown in Figure 3. 

Application of the box method to G 0 ( ip ), the duct face contribution to the far-field 
velocity potential, enables the integral equation to be written as 

M 

g„(0)=£ 

5=1 


. 


) J m (“Jr *nSin ^ 


(11) 


where the locations at the control points are 

r 5 = *3-1+— J - 2 - 3 — f °r 0^2 

and 


r 5 = 


At? 


+ a 


for j-i 


For the duct wall components, the equations become 


n r ex L ir( n '~ s> 

a,, 

5=1 ■Mss, . 

**-s0 


(12) 


where the locations of the control points are 


z. .= z, . ^l,j+ ^j-1 
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{13) 


G*(W=* 


sa 


c 



where the locations at the control points are 


<4z 


z 2,j = z 2,3-l+' 


2,3 +Az s,i-i 


Equations {11) through (13) are integrated numerically. The numerical technique used is 
an even order Gaussian quadrature scheme. The integration for each annular segment is 
performed by an iterative method. An initial even order quadrature is chosen (usually 
2nd or 4th order) and the integral evaluated. The order is increased by 2 and the 
integration repeated. 

The two values of the integral are compared. If they agree to within a specified error, the 
next segment is integrated. If agreement is not obtained, the order is again increased and 
the procedure repeated up to a specified maximum order. At this point the last 
calculated value of the integral is used and the integration procedure passes to the next 
segment. 

The quadrature method used is the Gaussian formula for arbitrary intervals. 



b 

f(x)dx 


b-a 

2 


NG 


0=1 


(14) 


This relation states that the integral may be approximated by a summation of IMG (the 
order of the quadrature) evaluations of the function at a prescibed abscissa, X:, which 
are multiplied by the associated weight, Wj. The abscissae and weights of the Gaussian 
quadrature may be found in Reference 3. The referenced abscissae have the normalized 
range -1 < y < 1 and are related to the arbitrary interval abscissae by 


x. 

J 


(b+a) (b-a) 


The numerical equations for the evaluation of the directivity functions are 


M NG 

0=1 *=1 


(15) 
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The quadrature points on the duct face are 


= r j + 


ir j x k 


N pi-f (»! 4-“) 

G, m = sb J m (f bSin <P )£) #,(*! •,) 2 u ’ V z l,j (1B) 

3=1 k=l ^ 


The quadrature points on the outer duct wall are 


Z l,3k = z l,j‘ 


K2 KG q^( U 2 i' S ) 

G 2 (^}= v '* z ” J (17) 


j=l k=l 


r"2 s j 


The quadrature points on the extended centerbody duct wall are 


Z 0 -Ur z o 


dZ 2,3 X k 


2, jk g 
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4.0 


RESULTS 


4.1 Application 

The numerical techniques developed in Section 3 are applied to the mathematical 
representation of the radiation problem developed in the Acoustic Source Distribution 
Program, Reference 1. Duct segmentation parameters and corresponding source 
distributions are input to this program to determine the far-field radiation characteristics 
of the duct configuration. Far-field directivity patterns are determined for each mode 
and then combined to yield the total noise signature. In the event both an inlet and fan 
duct are to be modeled, their respective far-field directivity patterns can be calculated 
and then combined to yield the total noise signature of the nacelle. 

An example duct radiation problem is presented to illustrate the output of this program. 
A source distribution which was generated by Reference 1 for a plane wave velocity 
distribution on the face of a circular duct {example problem 1) was input to this 
program. The resulting normalized far-field radiation pattern shown on a polar plot in 
Figure 4, for a circular duct with kb = 1.0, is compared with results obtained by a 
numerical calculation based on the closed form solution given by Reference 4. 
Comparison of the two results indicates that for this case the analysis does not predict 
the sideline portion of the directivity pattern very well. 


4.2 Limitations 

There are several limitations to this analytical program. First the effects of a mean flow 
within the duct on the radiation of sound has been neglected. Although application of 
this analysis to a nacelle radiation problem does not yield an exact model, it is 
anticipated that this simplification will not significantly affect results for noise 
propagation through a locally subsonic flow field of an inlet. However, radiation through 
an exhaust or jet flow will not be represented properly by the present analysis. 

The second limitation to the program is restriction of the specified velocity distribution 
on the face of the duct to a Cos (m Q ) type of angular dependence. This assumption 
does not affect the capability of the program to combine different m angular modes in 
the far-field by superposition of the source distributions. It does, however, restrict the 
Duct Termination Impedance Program, Reference 2, to the extent that for a velocity 
distribution whose angular dependence is specified by a combination of different m 
angular modes a separate analysis must be made for each set of angular modes. 


PRECEDING PAGE BLANK NOT EILK0S 
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5.0 


CONCLUSIONS AND RECOMENDATIONS 


An acoustic radiation analysis has been developed that predicts the angular distribution 
of acoustic energy in the far-field that is radiated from an unflanged duct. In spite of its 
many advantages, this acoustic radiation analysis is limited in its range of application. 
For this reason, this analysis has been relegated to the role of a first generation analysis 
in a commitment to develop an analytical package which will predict radiation 
characteristics of a turbofan jet engine. As discussed in Reference 1, the present version 
of the analysis reaches optimum efficiency when it is applied to low frequency noise 
radiating from an infinite duct. This and the other limitation due to a Cos {rc\6 ) type of 
dependence of the velocity distribution are ar. inherent property of the "box method" 
numerical technique. 

A second generation version of the program could avoid these limitations by utilizing a 
collocation technique of evaluating the integrals of Equations (8} through (10). The 
primary feature of this technique is that a selected analytical expression is used to 
represent the source distribution. The time saving feature of the collocation technique is 
that it requires fewer evaluations of the integrals. The collocation procedure would 
enable the analysis to handle any combination of angular modes, a result which would 
remove the angular modal coupling limitation for the Impedance Program, Reference 2, 
and eliminate the need for a large number of angular boxes. 

In addition to the collocation procedure, a second generation version of the analysis 
should contain a method of accounting for a mean flow within the duct. This method 
should account for the effects of velocity gradients on the distortion of the radiation 
field by the jet 


J8B8G0HNG PAGE BLANK NOT FUMED! 
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Computer Program Variables 


NOMENCLATURE 


A 

B 

C(l| 

C2(l) 

CHK 

D{l) 

FFR 

G{I) 

GR{I) 

INK 

K 

MW 

M 

N 

N2 

NG 

PHI{I) 

POP(I) 

R 

RO, Rl 
S 

TT 

Tl 

W 

WDS 

XANG 

ZZ 


Centerbody Radius 
Outer Duct Radius 

Dimensionless Size of Ith Box on Duct Outer Wall 
Dimensionless Size of Ith Box on Extended Centerbody Wall 
Convergence Criterion 

Dimensionless Size of Ith Box on Duct Exit Face 
Far-field radius 

Directivity function of Ith Polar Angle 
Normalized radiation pattern of Ith Polar Angle 
Configuration Parameter 

No. of Polar Angles at which Far-Field Data is required 
Angular Mode No. 

No. of Annular Boxes 

No. of Cylindrical Boxes on Duct Wall 

No. of Cylindrical Boxes on Extended Centerbody Walt 

Order of Gaussian Quadrature 

Source Strength Distribution of Ith Box 

Polar Angle array in degrees 

Radial Coordinate 

Distance between Point on Duct and Field Point 
Sound Speed in Air 
Angular Coordinate 

Angular Coordinate Integration Variable 
Frequency Radians/Second 
Wavenumber W/S 

Nondimensional Angular Size for Control Box 
Axial Coordinate 
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The basic operation of the Directivity index Computer Program is described and the 
input/output formats are listed. A flow chart of the program is presented in Appendix 2 
and listings of the main program and its subroutines are presented in Appendices 3 thru 
5. The input preparation and output format of the program used for the example 
problem in Section 4.1 is presented in Appendix 6. 

iNPUT 

The input data for the program consists of 16 variables whose symbols and card formats 
are described below: 

CARD OR 

CARD SET FORMAT DATA 

1 3612 K Number of angular points, , at which the 

directivity index will be calculated 

M Number of annular rings used to represent the 

face of the duct 

N Number of cylindrical rings used to represent 

the outer duct wail 

N2 Number of cylindrical rings used to represent 

the extended inner duct wail (if present) 

INK Indicator for extending the inner duct wall past 

the duct face. O does not extend the wall, 1 
extends it to infinity. 

ND Maximum number of passes to be made 

through the quadrature iteration loop. ND is 
related to the highest quadrature order NG by 
ND = NG-2 
2 

MW Angular mode number 

2 6E12.6 A Inner duct radius 

B Outer duct radius 

W Circular frequency, radians/second 

S Speed of sound in units consistent with those 

of A and Ei 

CHK Convergence criteria for numerical integration 
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CARD OR 
CARD SET 

FORMAT 


DATA 



FFR 

Separation Distance 

3 

6E12.6 
1 = 1, M 

D(I) 

Annular ring box widths. These widths a.e 
normalized by the outer duct radius, numbered 
D(1) at the duct centerline {or adjacent to the 
centerbody wall-duct exit plane intersection/ 
and 1 increasing to D(M) at the lip of the duct. 

4 

6E12.6 
1 = 1, N 

C(l), 

Outer duct wall box lengths. These lengths are 
normalized by the outer duct radius, numbered 
with C(i) closest to the duct exit plane and 1 
increasing in the negative z direction. 

5* 

6E12.6 
1 = 1,N2 

C2{l), 

Extended centerbody duct wall box lengths. 
These lengths are normalized by the outer duct 
radius, numbered with C2(l) closest to the duct 
exit plane and t increasing in the positive z 
direction. 

6 

6E12.6 
1 = 1, M 

PHIO(I), 

Source distribution on duct face. Method of 
numbering this source distribution is analogous 
toD(l). 

7 

6E12.6 
1 = 1, N 

PHll(l), 

Source distribution on outer duct wall. Method 
of numbering this source distribution is 
analogous to C{1). 

8* 

6E12.6 
1 = 1,N2 

PH 12(1), 

Source distribution on extended centerbody 
wall. Method of numbering this source dis- 
tribution is analogous to C2(l). 


Card Sets 5 and 8 are needed only if I NK = 1 
PROGRAM FLOW 


Program operation consists of two basic loops, one of which calculates the directivity 
function and the other which determines the normalized radiation pattern. An initial 
value of the angle $ is calculated. Equations (14) through (16) are solved to yield the 
exit plane and duct wall contributions G-j, G 2 , G 3 to the far-fiefd directivity function. 
The maximum absolute value of the directivity function G is found and used to 
calculate the normalized radiation pattern as follows: 


GR « 20 log 


( 


® ) 

G mox ' 
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Calculation of the directivity function requires two subroutines, GAUSS and BESJY. 
GAUSS is the subroutine that provides the weights and abscissae for the Gaussian 
quadrature described in Equation {14). BESJY supplies the value of the Bessel function 
J m for the argument XARG. Listings of GAUSS and BESJY are presented in Appendices 
4 and 5, respectively. 


OUTPUT 

The output format consists of three groups of data: 

® Input Data — All pertinent input data are printed out. 

* Directivity Function - The directivity function array and its associated 
array of angles, \p , in degrees, are printed out in column format. 

• Radiation Pattern — The normalized radiation pattern array and its 
associated array of angles, (p , in degrees, are printed out in column format. 
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APPENDIX 2 

PROGRAM FLOWCHARTS 




BLANK** 0 * 
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C NEW DIRECTIVITY INDEX PROGRAM 

DIMENSION C(40), P0PC180) , 0(40) , C?(40) 

COMPLEX *08 

1 PHI0140), PH 1 1 (40 5 » PH 1 2 140 ) t GU80) 

COMMON /GEN/ W, St CHK, NN, ND 
ICO FORMAT! 3612 J 
200 FORMAT! 6F12.61 

101 FORMAT t*l ENTERING THE DIRECTIVITY FACTOR PROGRAM . */// ) 


$ £ # # 


* # * * * # # * * * 


WRITE 16t101) 
READ { 5 T 100 } 

READ { 5 » 2001 
READ ! 5,200) 

READ ( 5,200) 

READ 1 5,200) 

IF ( INK *GT 
2000 READ15, 200 

READ { 5,200} 

IF { INK .GT. 0) 
CALL UFDIR (K 
1 ,POP, N 2 , 


K , H, N, N? , INK, ND, 
A, B, W, S, CHK, FFR 
! C ! 1 5 , 1 -I ,N } 

ID! I) , I =1 ,M) 

(POP! I), I =1 ,K 1 
0) RF AD (5,200) 

END-9000) (PHI 0( I) , 1 = 1 ,M] 
(PHIK I ) , I = 1 , N) 

RF AD ( 5,200) 
N, INK, MW, PHIO 


MW , NN 

!C2(I) , I- 1 ,N2 ) 

( PHI? (11,1=1 ,N2 ) 

PH 1 1 , PH 1 2 , A , B , D, C 

GO TO 2000 


, 


»Mf 
FFR ) 


9000 CONTINUE 
STOP 
END 


, C2,G 
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SUBROUT IMF UFDIR 

1 (K f M » N, INK t MWfr PHI 0 1 PHI1» PHI2, AD , BD » D, C, C2,G 

2 , PD P , N2 t FFR i 

C DIRECTIVITY FACTOR - UFDIR ZDM915 06 - 06 

REAL *8 AD t BD 

C XCNCt WTCND), XAR(NO), GR(K), GC K ) 

DIMENSION Ctl >,X{20) , WT 1 2 0) * PC1P 1 1 1 ,GR 1 170) , D(H t C 2 C 1 > 
COMPLEX *08 

1 PHIO(l), PHlltl), PH I 2 1 L ) i XAR(201t 

2 Gif G?t G3,Xl,X2,X3 f Ga80),DMPLX 
COMMON /GFN/ W, S, CHK, NN* ND 
EQUIVALENCE ( A ,AE 1 , t B , BE J 

C 

330 FORMAT I'O ANNULAR RING BOX WIDTHS 1 / { 8 XI P4F 1 5* 4 ) 1 
400 FORMAT C»0 OUTER DUCT WALL BOX LE NG THS ' / ( 8X1 P4E 1 5. 4 1) 

BOO FORMAT no INNFR DUCT WALL BOX LENG TH S * / 1 8 XI P4E1 6. 4 1) 

6C0 FORMAT no SOURCE DISTRIBUTION ON. DUCT FACE.'/ { 8X1 P4E 16.4) ) 

700 FORMAT t'O SOURCE DISTRIBUTION ON.OUTER DUCT WALL. '/ (8 X 1P4F 15.4)1 

SCO FORMAT CO SOURCE DISTRIBUTION ON. INNER DUCT WALL. * / 18X IP4F15.4) l 

1000 FORMAT { 3612) 

1010 FORMATt 6E12.6) 

1020 FQRMATL * P LIN DEGREES) = *, F13.5,1CX'G =',2E13.5) 

1030 FORMAT {' ** INTEGRATION OfO NOT CCNVFPGF. DIFFERENCE = 

1 F10. 4, 1 CHK =*,F10.4, ' R =*,F10.4) 

1050 FORMAT! 1H , 10X , 1 NUMBF R OF PSI POINTS = * , I 5 ,5 X T 1 A NGU LAP MHDF * , 

1* NUMBER t 151 

1060 FORMAT { 1H ,10X, ‘NUMBER OF ANNULAR RING SEGMENTS = *, I5/11X 

1 , ‘NUMBER OF OUTER DUCT WALL SEGMENTS = *, I5/11X, ‘NUMBER 0F» 

2 i* INNER DUCT WALL SFGMFNTS = *, I5/11X, ‘HIGHFST ORDER GAUSS* 

3 t * I AN QUADPATIJRF =*,I5) 

1070 FORMATt 1H ,10X,»IS INNER OUCT PRESFNT 0=ND,1 = YES * T I 5 ) 

1080 FORMAT! 1H ,10X,*DUCT OR “ * , F 1 0. 4 , * F T 1 1 5 X ,« DUC T l» =*,F10.4, 

l'FT'f/f 5X, ‘FREQUENCY =*,F10.4,* RAD /SEC * T 5X , ' SPEFD OF SOUNO=* , 
2F10. 4, 2X, ‘FT/SEC 1 /, 5X, * INTEGRATION CONVFRGFNCE NO =*,F10,4) 

1100 FORMAT!* P LIN OFGRFFS) -*, F13.5,10X, 'SPL = *, F13.5) 

12C0 FORMAT (»1*,?5X,‘0 IRFCTIVTTY FUNCTIO N»/l 
1300 FORMAT ni*,?5Xt*R AD IATION PATTER N* / ) 

1400 FORMAT t'l* T 25X 'INPUT DAT A*/i 
C 

C # £ * if if if if if $ * 

5577 FOPMAT (»C £UFDIR$ * T 8G15.7/ ( 10X , BG1 5 . 7 ) ) 

15 KP ••= K 
20 CO NT IN UF 
AF = AD 
RE - BO 

WRITE (6,1400) 

WR ITE( 6, 1P50J K ,MW 
WR ITE( 6, 1060) M,N,N2, NO 
WR ITEC 6, 1070) INK 
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WR I T E { 6 , 108Q } B,A,W,$,CHK 
WR ITS < 6,330} l DU) , 1=1, M) 

wr ire < 6 , 400 ) < c ( n , i=i,ni 

I F ( INK ,GT. 0} WRITF 16,500) (C2 U ) , 1 = 1 ,N2 I 

WRITE (6,600) ( PHIOt I), 1=1 ,M) 

WRITE (6,700} tPHIim,I = l,N) 

IF (INK, GT • C) WRITE (6,800) t PH 12 ( I ) ,1=1 ,N2 ) 

C - - BASIC INTTI L EZATION 

WDS = W /' S 

BD2 = B / 2,0 

SOW = S / W 

WBDS = W * 8 / S 
WADS = W * A I S 
AINK = A * INK 
WRITE (6,1200) 

DO 31 L = 1 , K P 

P = POP ( L ) / 57,295779^130823 FO 
IFUP.GT.l. 569), AND. (P.LT. 1.575)) P= 1,567 
27 SIP = S TN( P 1 
COP = CQS( P ) 

POP(L)=P* 57, 29577 

C LOOP INITIALIZATION 

SI = SDW / COP 
AS I = A # S I 
BSI - 8 * SI 
XARGB = WROS * SIP 
XARGA = WADS * SIP 
WDSSIP = WDS * SIP 

COMPUTE CONTRIBUTION QF DUCT FACE 

G1 = (0*C,0.0) 

no ioo i= i , m 

IF(I-l) 121,121,122 

121 R = D( I )*BD2 + A INK 

GD TO 122 

122 IL = 1-1 

R = R + (D( n+0( IL1 i*RD2 

123 XI = (0.0, 0,0) 

BD2D = RD2 * 0(1) 

C NO IS THE HtGEST ORDER GAUSSIAN QUAD USED. 

DO 150 KK=1,N0 
ILK = KK 

TLL =NN+2*CKK-1) 

CALL GALSS ( X, WT , I LL ) 

XAR(KK) ^ (Q.OrC.O) 

YAR = 0.0 

DO 155 11 = 1, ILL 

RQ - R - BD2D *XUI) 

XARG = RQ* WDSSIR 

CALL RESJY(XARG, MW, 3F1, PUM, DUM, OUM, I E RR ) 
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155 YAR = Y AR +8 F1*R 0 * WT(II) 

X AR { KK 3 = CMPLK ( Y AR*R * Pdl *.5, 0 o fE0) 

IFI KK-2 3150, 156 t 157 

156 KL = KK~ 1 

ABT = A BS( C ABS(XARIKK)} - C AB 5 ( XAR (KL ) ) 3 
IFI ABT-CHK ) 160,160,150 

157 IF(KK-ND) 156,153,160 

153 WRITE ( 6, 1030 3 ABT * CHK, R 
150 CONTINUE 
160 CONTINUE 

XI = XI + XARULK) 

100 G1 = G1+X1*PHI0U ) 

COMPUTE CONTRIBUTION OF OUTER DUCT WALL 

G2 = {O.C,O.OJ 

DO 200 1 = 1 1 N 

IF( I“I ) 22 I* 22 1» 222 

221 11 = -C ( I)*BD2 
GO TO 223 

222 IL = 1-1 

11 = ZZ— <C t I J+C C I L > )*BD? 

223 X2=(O.C,C.O) 

□0 250 KK=1,ND 
ILK =KK 

ILL=4+2*(KK-1 ) 

CALL GAUSS( X T WT, ILL 3 
XAR(KK)={0.0, O.C) 

DO 255 11=1, ILL 

Z=ZZ-B*CU 3*Xt I 1 3 /2.C 

RD=( FFR **2+Z*?-2. 0* 7*FFR*CDP )**0* 5 

R1 = ( RO—FFR 3 #W/ S 

CS=COS( R13/R0 

SI=SIN( R13/R0 

255 XAR(KK) =XAR (KK 1+CMPLXC S, ST )*WT( II) 

XAR{ KK) =XAR(KK)*B*C ( 13*0. 5 
IFIKK-2) 250,256, 257 

256 KL=KK~1 

ART= ABS( CABS ( XAR (KK 3 3-C AB S ( XAR( K L 3 3 ) 

IF ( ABT-CHK 3 260,2 60 ,2 50 

257 IF(KK-NC) 256,253,260 
253 WR I TF{ 6 , 103Q 3 ABT,CHK ,2Z 
250 CONTINUE 

260 CD NT INUF 

XAR G=W*EW= SIP/S 

CALL 3ESJY { XAR G , MW, B F2 , DUM, DUM, DUM, I FRR 3 
X2 = X 2 + X AR ( ILK )*BF2*FFR*B 
200 G2=G2+X 2TPHIK I 3 

G(LJ = G1 + G2 
IF ( INK .NE. 1 ) 
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C 

c 

c 


COMPUTE CONTRIBUTION OF INNER DUCT WALL 


50 G3 - {0,0,0,05 
DO 300 K1»N2 
TF{ I-I ) 321, 321, 32? 

321 ZZ=C2{ I KB/2,0 
GO TO 323 

322 IL = 1-1 

ZZ=ZZ + IC2m+C2IIL) KB/2.0 

323 X3= ( 0,0,0,01 
on 350 KK-l , ND 
ILK=KK 

ILL =4+ 2* ( KK-l ) 

CALL GAUSSl X, WT, ILL I 
XARI KK) =( 0,G,0.0) 

DO 355 1 1 - 1 , ILL 

Z = 7 7 +FKC2I I)n(lTl/2.0 

RO=C FFR**2+Z*Z~2. 0* Z*?FR*CrP K* 0. 5 

R 1= ( RO— FFR KW / S 

CS=COS{ Rl) /RO 

SKSINt RD/RO 

35 5 XAR(KK ) =XAR t KK ) +CMPLX I C S, SI ) * WT U I ) 

XAR { KK ) =XAR { K K K8*C ?t I KO. 5 
IFIKK-2) 350,356,35? 

356 KL = K K— 1 

AB T = ABS (CAB S( XAR IKK)) -CABS! XARI KLH I 
I F ( ART-CHK ) 360,360, 350 

357 IF(KK-ND) 356,353,360 
353 WR I T El 6,1030) ABT,CHK,ZZ 
350 CONTINUE 

360 CONTINUE 

XARG=W*A*SIP/S 

CALL BESJY I XAR G, MW, BF3, DUH, DUM, OUM, IERR) 
X3 = X 3+X AP t ILKKBF3*FFR*A 
300 G3=G3+X3*PHI2( I ) 

GtL) = GIL) + G3 

CALCULATE NORMALIZED RADIATION PATTERN 

30 WRITE! 6,1020) PGP(L),G(L) 

31 CONTINUF 
J-l 

* GMAX = C ABSI Gt 1 ) ) 

GR I 1 KGMAX 
60 J = J + 1 

GR ( J ) = C ABSIGI J ) ) 

IF(J.EO.KP) GO TO 70 
IFC GR{ J ),GT,GMAX) GO TO 75 
GO TO 6 C 
75 GMAX=GR(J) 
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GO TO 6 C 

70 IE( GRU l.GT.GMAX) GMA X-GR ( J 5 
WRITE ( 6*1300} 

DO 90 I - 1 1 KP 

90 GR{ I J=ALQG10(GR llj/GMAX }#2T. 

WRITEt 6,1100) {POP (I )tGR (I) ,I=1*KP) 

RETURN 

END 
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APPENDIX 4 


GAUSS SUBROUTINE LISTING 


PAGE BLANK HLMEIJ 



SUBROUTINE GAUSS t X, WT , M) 

C STORE DP TEST IN 81 ZDH915 11 - BT 

CD IMPLICIT REAL*8 {A-H,D-Z1 
DIMENSION XI 121 ,WT 112} 

2FIM-4)10,U,12 

10 X(l)=-0. 577350269189626 

xi 2 j= -xm 

WTll)=l.G 
WT (21=1.0 
GO TO 20 

11 XI n=-0,861136311594053 

XI 2}=-0.33598lC43584856 
X{3)=- XI 2 1 

x 1 4 > =- x m 

WTU 1=0. 347854845137454 
WT( 2) =0-652 145154862546 
W T t 3 I -WT t 2 1 
WTI4)=W7I 1} 

GQ TO 20 

12 IFIM-8113, 14, 15 

13 X( 11=-C. 932469514203152 

XI 2 J =-0-661209386466265 
XI 3) =-0-2 3861 9186083 19 7 
X I 4 J =- X t 3 I 

X I 5 ) =— XI 21 

xi 6 )=- xm 

WT I 11 = 0-1713 24492379 170 
WTI2) =0-360761573048139 
WT ( 3) =0-46 79 139 34572 691 
WTI4)=WTI 31 
WT ( 5 1 =WT l 2 ) 

WT t 6 J =W T t 1) 

GQ TO 20 

14 Xt II =— 0-960 2898 5649 75 36 
XI 2 )=“0- 796666477413627 
XI 3 )=-0. 525532409916329 
X ( 4 }=— 0- 18343464 249565 0 

X I 5 } =- X I 4 1 

X I 6 1 =“ XI 3 } 

xm =- xi 2 j 

X I 8 J =— Xtu 

WT ( 1 ) = 0- 1012285 36290376 
WTI2 1=0.222381034453374 
WT I 3 1 = 0.313706 645 8 7788 7 
WTI 4 1=0.362683783378362 
WT I 5 1 =WT ( 4 1 
WT I 6 1 =WT ( 3 1 
WTI 7 1 =WT I 2 ) 

WT t 8 ) =WT{ 13 
GD TO 20 

15 IFIM-121 16, 17, 18 
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16 Xt l)=-C. 973906528517172 
Xt 2 J=-0.S650633666B8985 
X( 3 >=-0.679409568299024 
X( 4) =-0.4 33 395 39412924 7 
Xt 5) =-0.148B74 338981 631 
X ( 6 1 =- X ( 5 J 
X(7I=- X ( 4 1 

x t a ) =- x { 3 } 

X { 9 ] —— X t 2 1 

x t io3=- xm 
Wit 11=0.066671344308688 
WT (21=0.149451349150581 
WTO 1 = 0. 219086 362515982 
WTt 41 =0.269266719309996 
WTt 51 = 0.295 524 2247147 53 
WTt 61=WT(5| 

WTt 71=WTt41 
WT t 8 3 =WT { 3 ) 

W T 1 9 1 =WT (21 
WTt 10 ) =WT l 1 ) 

GO TO 20 

L 7 X{ 1 l=-0. 981560634246719 
Xt 2) =-0.904117256370475 
Xt 3)=-0. 769902674194305 
X( 4} =-0.5873 17954286617 
Xt 5 l=-0. 367831498918 180 
X t 6 ) =-0 ■ 1253334C851 1469 
X t 7 1 =- X 1 6 1 

X ( 8 1 =— X { 5 1 

X t 5 } =- Xt 4 1 

Xtl01=- Xt 33 

Xtlll=- XI 23 

X t L2 1 =— X(U 

WTtl) =0.047 175336386512 
WTt 21 = 0.106939325995318 
WTt 31=0.160078328543346 
WTt 4 1=0. 203 167426723066 
WTt 51 = 0.23349 25 36 538355 
wTt 61 =0.249 147045813403 
WTt 71=WTC6) 

WT { 83 — WTt 51 
WTt 9)=WTt41 
WTt 10)=WTt31 
W T 1 11 1 =W T t 2 1 
WTt 12)=WTtl) 

18 IF t H — 16 5 

19 XI 1)=-C. 986283808696812 
XI 2 3 =-0.92 8434883663 574 
Xt 3) =-0.827201315069765 
Xt 4 )=-C. 687292904 81 1685 


19, 21, 20 


GO TO 20 
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X£ 5) =-0.5 1524B636358 154 
X( 6 3=-0. 3 191123 6892 7890 
X£ 7 }=-0. 10805494 8707344 
X £ 8 J =- X ( 7 1 

X£5)=- X { 6 J 
X £ 10 3 = -x(5i 

X(ll)=- XUJ 

X { 121=— X £ 3 ) 

X t 13 3=- X 1 2 J 

X £ 1 4 ) — — X £ 1 ) 

WT £ 11 = 0.035119460331752 
WT £ 2 3=0.0801580 87155760 
WTt 3 3=0. 121 518 570687903 
WT£ 4) =0,1572031671 58194 
WT£ 5 ) = 0. 185538 39 747793 8 
WT£ 6 3 = 0.205158463721296 
WT£ 7 3=0.215263853463 158 
WT < 8)=WT{7) 

WT £ 9)=WT£6) 

WT ( 10 ) =WT C 5 3 
Wtl 113=WT£43 
WT £ 123=WT£ 3 } 

W T £ 13 1 =WTt 2 3 
wT £ 14 ) =W T £ 1 J 

21 X( 1 )=-0. 989400934991650 
X( 23 =-0.944575023073233 
X( 3)=-C.86 563 1202387832 
XI 4)=-0. 755404408355003 
Xt 5 3 =-0.617876244402 644 
X £ 6 3 =— 0 .A 580 1677765722 7 
XI 7 )=-0. 281603550779259 
Xt 8 3 =-0.0950 1250983763 7 
X { 9 3 =— X £ 8 3 

X £ 10 3= -X £ 73 

X£ 11)=- X( 6 J 

X t 121=- X ( 5 } 

X £ 13}=- X£ 

X t 14}=— X t 3 } 

X £ 15)=- X t 2 ) 

X £ 16 ) =— X(1J 

WT£ 13=0.027152459411754 
wT £23=0.062 25 3523938648 
WT £ 3 3=0.09515851 1682493 
WT{ 43=0. 124628971255534 
WT£ 53=0.149595988816577 
WT {63=0.1691565 19395003 
WT£ 7 3 =0.18260341504492 4 
WT £8 3=0.18945 0610455069 
WT £ 9 3 =WT £ 8 3 
WT £ 10 ) =WT£ 7 } 


GO TO 20 
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WT( 11) 

wr ( 12) 

WT 1 13 ) 
WTU4) 
art is > 
hT( 16} 
20 RETURN 
END 


=WTI 6 ) 
=WT t 5 ) 
= W T ( 4 ) 
S WT{ 3 ) 
=V» T t 2 > 
=* 7 ( 1 ) 



APPENDIX 5 

BESJY SUBROUTINE LISTING 


PRECEDING PAGE BLANK NOT FUMED 
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SUBROUTINE BES JY £ X,NO, B J,BY, BJ J, BYY* I ERR) 
REAL J(350),J0, Jl, JOO 
DIMENSION Y ( 60 } 

PI-3.141593 
IFIX.GT.O.) GO TO 90 
I ERR = 0 
GO TO 80 
90 CONTINUE 
I ERR= 1 
FAC = 1*5 

IF £ X *GT. 50.0) FAC = 0.5 

NX =3 • 9+FAC+X+NO 
J( 2*NX+2}=0. 

Jl 2*NX+1 )=l*E-29 
C=0. 

J0=0. 

J1=0. 

Y0 = 0 * 

Y 1 = 0 • 

CO 20 NB= 1 » NX 
N=NX+1-NB 
RN = N 

J £2*N)=2.*£2.*RN+1 }*J£2*N+1)/X-J (2*N*2) 

J ( 2*N- l ) =4 ■ *RN$J t 2*N i /X- J t 2*N+ 1 ) 
Y0=Y0*l-l)**N*J{2*N) /RN 

Yl=Yl+£— 1 £ 2^ N+ 1. )*J£2*N+1) /RN/(L.+RN) 
20 C=C+2.*J£2*N> 

JOO=2.*J ( U/X-J£2i 
C=J00*C 

Y0 = I ALOGt X/2. ) + . 57721 566 )*J00-2.*Y0 
Yl = l ALOG£ X/2. )-. 42278434 ) *J £1) -J OO/X-Y 1 
Y0=.63661<5 8*Y0/C 
Yi=.6366t98*Yl/C 
J 0 = JOO/C 
Y I l ) =Y 1 

Y£ 23-2. *Y£ 1 J/X-YO 

NNX=N0+2 

DO 30 1=1, NNX 

Y £ I + 2)=2.*FL0AT I l + l ) *Y £ I +1 1 /X- Y£ I J 
3C JUNIII/C 

IFtNO.EQ.OJ GOT 0 70 
BJ=J( NO) 

BY=Y{ NO) 

NP=N0+1 
BJ J =J { NP 3 
BYY=Y£ NP 1 
GQ TO 80 
70 B J = JO 
BY = YO 
BJ J=J ( 1) 

BYY=YUJ 
80 CONTINUE 
RETURN 
END 
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APPENDIX 6 


EXAMPLE PROBLEM 



891 6400Q00080C 230 8 

1,0 343*0 343,0 .0001 100, 

2* 50000D- 02 2*500000-02 5.00000D-02 5*000000-02 5.00000D-02 5.00000D-0? 
7*500000-02 7# 5CCCQD-02 7*500000-02 7.50000D-02 7.50000D-D2 7.50000D-02 
l.OCOOOn-Ol 1*000000-01 1*000000-01 1 * 000000-01 1*0G000D~01 1*00000D-01 
l.COOOOD-Ol 1 *000000— 01 1*000000-01 1*250000-01 1*250000-01 1.250000-01 

1.250000-01 I. 250000-01 1.250000-01 1. 250000-01 1.250000-01 1.250000-01 

1.250000-01 1.250000-01 1*250000-01 1.25000D-01 1.250000-01 1*500000-01 
1*500000-01 1.500000-01 1.500000-01 1.500000-01 

1.250000-01 1.250000-01 1.000000-01 7.50000D-02 5.00000D-02 2.500000-02 
2.500000-02 5.COCOOD— 02 5.00000D-02 7.500000-02 7.500000-02 7.500000-0? 
5. Q0GO0D-02 5.CC0CCD-02 2.500000-02 2.500000-02 


2*0 

4.C 

6.0 

8 * 0 

10.0 

12.0 

14*0 

16. C 

18.0 

20.0 

22.0 

24.0 

26.0 

28.0 

30* 0 

32.0 

34.0 

36.0 

38. 0 

40.0 

42. 0 

44,0 

46.0 

48.0 

50.0 

52. C 

54.0 

56. 0 

58,0 

60.0 

62.0 

64. C 

66 . 0 

68 . 0 

70.0 

72.0 

74.0 

76.0 

78.0 

80*0 

82*0 

84.0 

86*0 

88.0 

90. 0 

92. 0 

94.0 

96.0 

98.0 

1CC.C 

102.0 

104.0 

106.0 

108.0 

no.o 

112. C 

114. 0 

116.0 

118.0 

120.0 

122.0 

124. C 

126.0 

128.0 

130.0 

132.0 

134*0 

136. C 

138. 0 

140*0 

142.0 

144.0 

146.0 

148.0 

1 50.0 

152.0 

154.0 

156.0 

158.0 

16C.0 

162.0 

164. 0 

166.0 

168.0 

170.0 

172. C 

174. 0 

176.0 

173.0 



-6.051 150-0? 4 .396230-02-7.070610-02 4.862960-02-9.074830-02 5.779730-02 
-1.164650-01 6.955110-02-1.427180-01 8. 1 54 490- n 2~l. 6277 80-01 9.070B10-0? 
-1,734120-01 9 .734920-02-2 .C60 10D-01 1 . 1 04570-01 -2 . 5 2 7 1 2 0-0 1 1.317960-0] 
-3 . 338090-01 1 .68B75D-01-4. 855940-01 2.383490-01-7.459170-01 3.576790-01 
-1.101920 OO 6.209690-01-1.360760 00 6. 3 85 P3D-0I -1 . 6846 PD-0 1 8.127090-0? 
9,232 J2D 00-4.292080 00 

-5,246290 01 2.439520 01 3.305650 00-1.525480 00 9.598000 00-4.440850 00 
6.309650 00-2.905790 on 3.51748D 00-1.605960 00 1.942270 00-8.739000-01 
9.676640-0.1-4.221780-01 4.561 170-01-1. 866000-01 2 . 3 3 5 75 D-01-8 . 5 3 307P-02 
1.247160-01-3.677370-02 6.554990-02-1.116980-02 3.098380-02 3.159380-03 
7.029190-03 1 .244690-02-8.372440-03 1.768620-02-1.633480-02 1.967140-02 
-2.006600-02 1.987230-02-2.134460-02 1 . 908?on-02 -2 , 1 l 9720-02 1.774860-0? 
-2 . n 24 560-0? 1 ,61 5220-02-1 .890140-92 1.446 72 0-02-1.7409 2 0-0? 1.280810-02 
-1.579470-02 1. 101020-02-1. 39293D-9? 9. 1 59 830-03 -l . 2 36 1 1 0-T 2 7.506910-03 
-1.101260-0? 6.047700-03-9.867730-03 4.762740-03-8.895180-03 3. 630740-03 
-8.061170-03 ? .64105D-03-7. 338820-03 1.770590-03-6.691520-03 1.004170-03 
-6 . 107440-03 3.311 800- 04-5 . 560 72D- 03-? . 6 87300-04 -5 . 03 78 40-0 3-7 . 98 5 8 1 0-04 
-4.531 75 0-O3-1 .27 1280- 03-4 . 0 3 1 0 ID- 03- l . 7 Cl 690-03-3.46 8190-0 3-2. 148 880-0 3 
-2.82 5 96 0-0 3-2 .64320 0-03-2 , 1 30 0 30-03-3. ? 306 80-93 -l . 2 7 39 80-0 3-4 . 2683 1 0-0 3 
3.54277 D- 04- 8. 183340-03 
FNC 
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[ N P U T 


DATA 


NUMBER DF PS I POINTS = 89 ANGULAR MODE NUMBER * 0 

NUMBER OF ANNULAR RING SEGMENTS = 16 

NUMBER OF OUTER DUCT WALL SFGMFNTS = AC 

NUMBER OF INNER DUCT WALL SFGMFNTS = C 


HIGHEST OKDFR 

GAUSSIAN QUADRATURF = 8 


IS INNFP DIJCT 

PRESENT 0=N0 

♦ 1 = YF S 0 


DUCT OR = 

l-CCOOFT DUCT I R - 0.0 

FT 

FREQUENCY = 343.0000 RAD/SFC 

SPEED OF SOUND- 

343.0000 FT 

INTEGRATION CONVERGENCE NO = 0, 

0001 


ANNULAR RING BOX WIDTHS 



I . 2500E-01 

I . 250 OF -01 

1. 0000F-C1 

7*5 000E-02 

5.0OO0F-02 

? .500CF-02 

2. 50COF-C2 

5.0000E-02 

5 • CO COE- 02 

7 . 5000F-02 

7. 50 COE- C? 

7 .5000 E— 02 

5.P000E-02 

5 ,r 00CF-02 

?. 50 OOF— C2 

2.5 00OE-02 

GUT FR DUCT WALL BOX LENGTHS 



2. 50 OOF-02 

2 - 5OO0E-02 

5. OPCOE-C2 

5.0000E-02 

5.COOOE-02 

5.0000F-02 

7. 5rrop-C2 

7.5000F-0? 

7.50 00 E~0? 

7.50H0E-02 

7. 500OE-02 

7.5 TO0E-02 

l.ooooE-m 

l.f eocF-oi 

1. OOOOF-Ol 

l .OOOOF-Ol 

1.00 00F-01 

1 .nnooE-01 

1. O0O0E-CI 

i.OOOOE-Ol 

1 .onooE -01 

1 .250CF-01 

1, 25C0F-01 

1 .2500E-01 

1 . 250nE-OI 

1 .250GF-01 

1.25 OOF- Cl 

1 .2500E-01 

1*25 00 F-0 1 

1.2500E-0I 

L. 2500E-O1 

1.2500E-01 

1 .25D0E-01 

I .2500F-0I 

1. 25C0F-C1 

1 .5HOOF-01 

1. 5000 F-01 

1.5000E-01 

1. 50 OOF- 01 

1.50COF-01 

SOURCE UISTR I BUT I UN ON 

.DUCT FACE* 



-6.0511E-02 

4 .3963F-Q2 

- 7.P706F-C2 

4.8630F-02 

-9.07A9E-02 

5*770 7E— 02 

-1. 1646E-C1 

6.9551E-02 

-1 .4272E-0I 

8.1545F-02 

-1. 6278F-C1 

9 .0 708F-02 

-I. 7BA1F-01 

9.7 849E-02 

-2. 06P1E-C1 

1. 1046E-O1 

-2*52 7 lt-01 

1 .3 180F-0 l 

-3. 33 8 IF- Cl 

1.6387F-01 

-A. 85 59 E- 01 

2.3835E-01 

-7. 4592F-C1 

3.57A8E-01 

-L. 1019F 00 

5.2C97E-01 

-1. 3607F CO 

b . 3858E-01 

-1.6R46E-0I 

8*127 IF— 0? 

9. 2323E CO 

-4.2921E 00 
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SOURCE DISTRIBUTION UN. OUTER DUCT WALL. 


-5. 246 3 F Oi 

2 .4395E 01 

3.3056F rn 

-1.5255E 00 

9.5980F no 

-4.4408F 00 

6 . 2096E 0^ 

-2.9058F 00 

3 - *5 1 75T on 

— 1 .606 CE 00 

1.9423F 00 

-8.7390F-01 

9.6756 E— 01 

—4 .221 BE— 0 1 

4. 5612E-C1 

— 1.8 66PE-0 1 

2 . 33 57E-01 

— 8*5 33 IE — 02 

1 • ?4 72 £- Cl 

- 3 .6 7 74 E-02 

6.5550 C-02 

-1*117 CE-02 

3.C984F-C? 

3.1 594F-03 

7.0292F -Oi 

1 • 2 44 7F— 02 

- 8 . 3724C-C3 

1 .7686F-02 

-1. 63356-02 

1 .967 IF— 02 

-2. 0066E-C2 

1 *98725—02 

-2. 1344F-02 

1 .9C82F— 02 

-2. 1197E-C2 

l .7749E-02 

-2.0246F-0? 

1 .6152F-0? 

— 1.P901F-C2 

1 .446 7F-02 

- 1 • 74 09 E-02 

1 .2 808C-02 

-1. 5705E-C2 

1 . 10 10 E -02 

-1.3929E-02 

9.1598E-03 

- 1. 2361E-02 

7 . 5 06 9 E— 03 

-1 . 1013 E— n? 

6.O477E-03 

- 9 . 86 77C-C3 

A.7627E-03 

-8 . 8952E-P3 

3.6307E-03 

-8. 0612F-C3 

2.6410F-03 

— 7 . 33 88 E— 03 

1 • 7 70 6F-03 

- 6 . 691 5E-03 

I .0O42E-03 

-6.10 74 F - 03 

3.311 8F-04 

-5. 5607E-C3 

-2.6H73E-04 

-5.O37RE-03 

-7.9858E-04 

-4. 631 7F-C3 

-1 .2713E-03 

—A* 03 10F-0i 

- 1 . 7017C-03 

-3. 46H2E-C3 

-2.1 489F-03 

8 260 F- 03 

-2.6432F-03 

13P0E-03 

-3.2307E-03 

-1. 2740F-O3 

-4.2 68 3F-03 

3. 542HE-C4 

-H. 1 883F-03 
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D I R E c 

T I V I 

: r y 

F 

UNCTION 



p 

{ in 

DEGRFFS ) 

- 

0.200 COE 

01 

G 

= 

0. 1 C775E-01 

0.18837E 

00 

p 

( IN 

OEGREFS) 


0.4 10 POE 

01 

G 

- 

0. 1C591E-01 

0. 18792E 

CO 

p 

{ IN 

DFGRFFS) 

- 

0.600C0F 

01 

G 

ss 

0,1 P30/F-01 

0. 18721 E 

00 

p 

( IN 

DFGRFFS } 

= 

0.800C0F 

01 

G 

= 

0.9E903E-02 

0. 1B619E 

no 

p 

{ IN 

DEGREES ) 

= 

0.1 00 OOF 

02 

G 

= 

0. 937 89 E-Q2 

0. 18493E 

00 

p 

( IN 

DEGREES } 


0 . 1 20 00 E 

02 

G 

- 

0.875C9F-0? 

0. 18339E 

00 

p 

{ IN 

DEGREES ) 

= 

0.140COF 

02 

G 

- 

0. 80006 F-02 

0. 1R156E 

00 

p 

t IN 

DEGREES ) 

- 

0 . 1 6000 F 

02 

G 

tt 

0. 7 IS 55 F— 02 

0. 17951 F 

00 

p 

( IN 

DFGREFS ) 

= 

O.IBOCOE 

02 

G 

St 

0.6 2056 F— 02 

0.17720E 

00 

p 

I IN 

DFGRFFS } 

= 

0.2 0000 E 

02 

G 

:= 

0.51766E-02 

r . 17472E 

00 

p 

( IN 

0 FOR EES ) 

= 

0.?2f COE 

02 

G 

= 

0.40184E-02 

0. 17193E 

00 

p 

{ IN 

DEGREES 1 

= 

0 .240 COE 

02 

G 

- 

0.278UF-0? 

0. 16898E 

oo 

p 

( IN 

DEGREES} 

= 

0.2BQC0F 

02 

r 

= 

0. I 43 80 F-02 

0. 16584E 

00 

p 

l IN 

DFGREFS ) 

St 

0.2 ROOOF 

02 

G 

= 

0.56949E-04 

0. 162 56 F 

on 

p 

t IN 

DEGREES} 

= 

0 . 30C0OF 

02 

G 

St 

-0. 1 4178E-02 

0. 159 10E 

00 

p 

( IN 

DEGREES ) 

= 

0.3 20 00 E 

02 

G 

- 

-0. 2 5/04 E-02 

0. 1555CF 

00 

p 

{ IN 

QEGRCC 5 ) 

= 

0 .34000E 

02 

G 

= 

-0.4 6084 F— 02 

0. 15174E 

00 

p 

t IN 

DFCRFFS ) 

- 

0.360C0F 

02 

G 

= 

“0 . 6 3234-E— 0 2 

0. 14784 F 

00 

p 

{ IN 

DEGREES) 

= 

0.3 80 DOE 

02 

G 


-0.8C513E-02 

0. 1 4395 F 

00 

p 

I IN 

DEGREES } 


0 .40OC0E 

02 

G 

= 

— 0.9 8688 E— 0 2 

0. 13990F 

00 

p 

I IN 

DEGREES ) 

= 

O.420COE 

02 

G 

= 

-0.1 1724F-01 

C. I35 3 3F 

no 

p 

( IN 

DEGRPFS ) 

- 

0.440COE 

02 

G 

= 

-0.1 3636E-0 I 

0. 1316HE 

no 

p 

l IN 

DEGREES ) 

= 

0 .46'JCOE 

02 

G 

St 

-0.1 6568F-0 1 

0. 12754E 

00 

p 

( IN 

DEGREES ) 

= 

A.4H000E 

02 

G 

- 

-0. 1 7549F-01 

0. 12336F 

00 

p 

( IN 

DEGRFFS) 

= 

0 . 5QOCOF 

0? 

G 

tr 
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